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Abstract. We explore a chameleon type of interacting dark matter-dark energy scenario in 
which a scalar field adiabatically traces the minimum of an effective potential sourced by 
the dark matter density. We discuss extensively the effect of this coupling on cosmological 
observables, especially the parameter degeneracies expected to arise between the model pa¬ 
rameters and other cosmological parameters, and then test the model against observations 
of the cosmic microwave background (CMB) anisotropies and other cosmological probes. We 
find that the chameleon parameters a and /?, which determine respectively the slope of the 
scalar field potential and the dark matter-dark energy coupling strength, can be constrained 
to a < 0.17 and (3 < 0.19 using CMB data alone. The latter parameter in particular is con¬ 
strained only by the late Integrated Sachs-Wolfe effect. Adding measurements of the local 
Hubble expansion rate Hq tightens the bound on a by a factor of two, although this apparent 
improvement is arguably an artefact of the tension between the local measurement and the 
Ho value inferred from Planck data in the minimal ACDM model. The same argument also 
precludes chameleon models from mimicking a dark radiation component, despite a passing 
similarity between the two scenarios in that they both delay the epoch of matter-radiation 
equality. Based on the derived parameter constraints, we discuss possible signatures of the 
model for ongoing and future large-scale structure surveys. 
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1 Introduction 

There is a growing body of evidence that the universe is currently undergoing a phase of 
accelerated expansion (e.g., [1-6]). This apparent acceleration is usually attributed to a dark 
energy with an equation of state w ~ — 1. The fundamental nature of dark energy, however, 
remains unknown. At present, a cosmological constant or vacuum energy with w = — 1 
appears to be simplest solution that can account for all available data. But dynamical scalar 
field models of dark energy [7-9] or f(R) modified gravity models [10-12] remain viable, 
albeit highly constrained, possibilities. 

Yet another interesting possibility are the interacting quintessence models [13]. Here, 
a scalar field responsible for driving the accelerated expansion interacts with dark matter 
and/or ordinary matter. 1 In fact, in the absence of any symmetry forbidding the interaction, 
it is quite natural to expect such couplings [17, 18]. An immediate consequence, however, is 
that the matter fields will perceive an additional long-range force on super-Mpc scales; in the 
case the scalar field couples to baryonic matter, any such “fifth-force” effect that may have 
trickled down to the sub-Mpc scales is automatically subject to solar system constraints [19, 
20] as well as bounds from laboratory test of gravity [21, 22], 

One possible way to avoid these local constraints is the chameleon mechanism. The 
basic idea is that because the scalar field interacts with matter, its mass is a function of the 
local matter density. In those chameleon dark energy models wherein the scalar field couples 
universally to all forms of matter, the high baryonic matter density on solar system and 

Ht has been shown that any f(R) model can be mapped to an interacting quintessence model by a 
suitable conformal transformation of the metric to the Einstein frame. In a sense, both f(R) and interacting 
quintessence are scalar-tensor modifications of gravity that introduce a scalar fifth force [14-16]. 
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laboratory scales causes the scalar field to become heavy, thereby hiding the fifth force from 
local tests. This screening mechanism switches off on matter-thin super-Galactic scales, and 
ideally the fifth force should invert its behaviour so as to drive the accelerated expansion of 
the universe [23-27]. However, a recent work [28] has shown that the range of the chameleon 
force in universally-coupled models is necessarily restricted to below the ~ 1 Mpc scale, so 
that on its own the chameleon field is in fact unable to account for the observed accelerated 
expansion. This constraint can nonetheless be circumvented if the scalar field interacts only 
with non-baryonic matter, and the chameleon mechanism as discussed above remains viable. 

In this work, we consider one such interacting dark matter-dark energy (DM-DE) model 
proposed in [29], in which the scalar field adiabatically tracks the minimum of the effective 
potential starting from deep in the radiation-dominated era up to the present epoch. This 
model is very simple in that the coupling between the dark matter spinor t/j and the scalar 
field <j) is of the Yukawa type, f (fy'i/j'ip, and the coupling function is chosen to be a positive 
exponential suppressed by the Planck scale, i.e., /(</>) = which is very common in the 

Einstein frame and emerges from many string theory models of compactification [30] ; the self¬ 
interaction potential for </> takes the runaway form V(4>) ~ </>“". Previously, a naive estimate 
has put an upper limit on a, a < 0.2, while /3, which controls the strength of the fifth force, 
was found to be practically unconstrained [29]. In the present work, we wish to confront 
this model with the most recent cosmological data, especially measurements of the cosmic 
microwave background (CMB) temperature anisotropies by the ESA Planck mission [5], and 
reexamine the model’s viability. 

Our second motivation comes from recent hints of a possible excess of radiation energy 
density during the CMB decoupling epoch from measurements of the CMB damping tail by 
the Atacama Cosmology Telescope (ACT) [31] and the South Pole Telescope (SPT) [32], as 
well as from the combined analysis of Planck data and the Hubble parameter inferred in 
the local neighbourhood [5, 6]. Conventional explanations include a thermalised populatioin 
of eV-mass sterile neutrinos (e.g., [33, 34]) or other light particles (e.g., [35]), relativistic 
decay products of heavy particles (e.g., [36-38]), as well as early dark energy [39]. However, 
ultimately, what the CMB anisotropies probe are the effects of the energy content on the 
evolution of the photon perturbations around the decoupling epoch [40, 41], To this end, 
chameleon models, in which the dark matter-dark energy interaction endows the dark matter 
component with a nonstandard time evolution that tends to delay the epoch of matter- 
radiation equality, may very well mimic the phenomenology of a radiation excess. 

The plan of the paper is as follows. We begin in section 2 with a brief review of the 
chameleon model, and present the relevant equations of motion. In section 3 we discuss the 
effects of the chameleon model on cosmological observables, and the parameter degeneracies 
expected to arise between the model parameters a and /3 and other cosmological parameters. 
We perform a fit of the model to cosmological data in section 4, and discuss the implications 
for the chameleon model parameters as well as the possibility of chameleon dark energy 
mimicking dark radiation. Using the results from section 4, we identify in section 5 possible 
observational signatures for future cosmological probes. We state our conclusions in section 6. 

2 Chameleon model 

We briefly review the chameleon model of [29] in this section, and present the relevant 
equations that determine the background evolution, the initial conditions, and as well as 
the evolution of the dark energy and dark matter perturbations. These equations will later 
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be embedded into the Boltzmann code Camb [42], in order to analyse the impact of the 
chameleon mechanism on cosmological observables such as the CMB anisotropies and the 
large-scale structure distribution. 


2.1 Background evolution 

Following [25, 43], we take the Einstein frame as the physical frame in which dark energy 
interacts with dark matter fields through a 0-dependent conformally rescaled Jordan frame 
metric gfw , whose generic form is given by 


g 


(0 

/.W 


- p 2/3;</)/Mpi 


( 2 . 1 ) 


where /% are dimensionless coupling constants, and M P1 is the reduced Planck mass. The 
matter spinor fields i/)® are understood to follow the geodesics of g\ll , while the original 
9fiu = diag (—1, a 2 , a 2 , a 2 ) is the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 
in the Einstein frame with scale factor a. The action takes the form 


S = 


/ 


d A x 


yf~~9 


Mr _ M! _ V W 


+ 


a $: 


( 2 . 2 ) 


where 1Z is the Ricci scalar, C m the Lagrangian density for the matter fields, and we have 
also added the kinetic and potential terms for the scalar field <j>. Varying the action with 
respect to (j), we obtain an equation of motion for the scalar field: 


0 + 3H<j) 


-vM + E 

i 


Jh__ 

Mpi 


Af)i<t>/Mp\ n ^vrp(i) 
e %) 1 \xv 


(2.3) 


where an overhead denotes a derivative with respect to the cosmic time t, H = a/a is the 
Hubble expansion rate, = (~VV — g^>)5 C m / 5g^ the Jordan frame energy-momentum 
tensor of the ith matter field, and we have assumed that the scalar field is homogeneous and 
isotropic. Limiting the coupling to the dark sector, equation (2.3) reduces to 

0 + 3 Hj> = -V 0 (0) - A e W/M P1 - dm 9 (2.4) 

IVlpi 


where p DM = s is the Jordan frame mean dark matter energy density. 


9 t lV -b (DM) 

To determine the evolution of /5 DM , conservation of energy—momentum = 0 in 

the Jordan frame gives 

Pdm + 3 [ H + —— <j )) p DM = 0 


M P 


(2.5) 


which has the solution 


Pom — 


~( 0 ) 

Pdm 




(0) 

Pom 


/3(34>+<f>o)/M Pl 


( 2 . 6 ) 


where 0o = 0(a = 1), and /3 d°m = Pom^^°^ Mfi is identified as the present-day physical 
dark matter density in the Einstein frame. Substituting equation (2.6) into (2.4), we see 
immediately that the scalar field dynamics in this model is controlled by an effective potential 


Kfi(0) = R(0) + 


(o) 

Pom 


e P((j)-rj>o)/Mpi 


(2.7) 
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where, for this work, we adopt a simple power law runaway potential, 

V(t) = M}(^y , ( 2 . 8 ) 

with a a positive constant that characterises the scalar field’s self-interaction, and the 
scalar held mass. This form of potential is common in many string compactihcation mod¬ 
els [30], and has the property that the attractor solution possesses a growing DM-DE cou¬ 
pling. In the limit a = 0, the model reduces to the case of a cosmological constant. 

The Friedmann equation for the full system, including other standard forms of energy 
densities, reads 


H\t)= ( =-L^ 

W [ a) 3 Ml 


P? + P^ 


a 


a* 


+ 


(o) 

Pb 


+ 


(o) 

Pom 




T ~ h V ((f)) 


(2.9) 


where p^\ p ® and p^ denote the present-day energy densities of photons, massless neu¬ 
trinos and baryons, respectively. Again, for a vanishing coupling (5 and a constant poten¬ 
tial V((/>), equation (2.9) reduces to the standard ACDM equation. Observe that the DM-DE 
interaction induces a nonstandard time evolution for the dark matter density, i.e., 

(o) 

p DM = _ (2.10) 

a 6 

Furthermore, while the equation of state for the scalar held alone is the canonical one, namely, 


Pp _ 4 > 2 /2 - V(<f>) 
Pt </> 2 /2 + V((f)) ’ 


( 2 . 11 ) 


where p$ and are, respectively, the scalar held’s energy density and pressure, taking into 
account the nonstandard dark matter density (2.10), the effective dark energy density (as far 
as the background expansion rate H(t) is concerned) is in fact 


(o) 

Poe = P* + ^ (eW-*° )/Mpi - l) = (1 - x) P<j> , 
which has an effective equation of state 

w eS = w# x (1 - x) _1 . 


( 2 . 12 ) 


(2.13) 


Where the coupling increases with time, such as in our case, x ^ 0 holds, thereby driving 
w eff to a value below — 1, and hence leading the dark energy to exhibit an apparent phantom 
behaviour. 


2.2 Attractor solution 


The attractor solution for chameleon models and its stability have been studied previously 
in [25, 43, 44]. It is set by the minimum of the effective potential, i.e., V^((f) min ) = 0, and 
adiabatic evolution simply stipulates that cf{t) traces the minimum at all times. For the 
potential (2.8), this means 



e PW)~<t>o )/Mpi 
a 3 {t) 


0 ! 


(2.14) 
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Figure 1 . Left: The effective potential V e g (solid black) for the model parameters (a,P,uj OM ,h) = 
(0.2,0.1,0.13,0.7). It is a composite of the original potential V(p) (dashed blue) and the coupling 
with dark matter (dotted red). Right: The corresponding p field as a function of the scale factor a, 
obtained from numerically solving equation (2.4) (dotted red), and from the analytical solution (2.14) 
(solid black). 


which can be solved using a simple bisection algorithm once the scalar field’s present-day 
value p o has been determined. The latter can be achieved by demanding that the field’s 
present-day potential energy dominates over its kinetic energy, so that ~ V{pp). Then 
from Kff (</>min) = 0 we find 


Po ~ _ « / h 2 _ \ 

P pW P V w dm / 


(2.15) 


where the last equality in terms of the reduced density parameter o; DM and Hubble parame¬ 
ter h follows from the assumption of a flat spatial geometry. An estimate for the mass scale 
of p field, 






1/4 


(2.16) 


also follows immediately. Thus, once w DM an h have been specified, it is possible to uniquely 
parameterise an adiabatically evolving chameleon model with only two parameters a and P 
in a spatially flat universe. 

Figure 1 shows the effective potential V eS (p) for choice of parameters (a, P,u UM ,h) = 
(0.2,0.1,0.13,0.7), and the corresponding evolution of p —the numerical solution of equa¬ 
tion (2.4) as well as the analytical solution (2.14)— as a function of the scale factor a. Note 
that we have chosen the initial value of the scalar field p m to be in the undershooting regime, 
i.e., 0 < 0 in < p™ in . This ensures the satisfaction of the slow-roll conditions, 


V l <£ » H, P 2 « 2 V(P), (2.17) 

at all times, and consequently the stability of the linear density perturbations [44]. 

An immediate phenomenological consequence of the attractor solution is that the effec¬ 
tive dark energy equation of state w eS as defined in equation (2.13), and consequently the 
cosmological background evolution, can have no dependence on the coupling parameter p. 
This is easy to see by first defining p = Pp/M Pl . Then the present-day field value is given 
by po ~ a(h 2 /u DM — 1) according to equation (2.15). Substituting this into the attractor 
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solution (2.14), we see straight away that <f>{t) does not depend on j5. It follows then that w eB , 
whose modification due to the DM- DE coupling comes in through the exponential e ^~^°, is 
also independent of /3. To measure /3 we need to consider its effects on the evolution of the 
inhomogeneities. 


2.3 Evolution of inhomogeneities 

The presence of a DM-DE coupling likewise modifies the evolution of the inhomogeneities in 
the universe. Working in the synchronous gauge whose line element is 

ds 2 = — dt 2 + a 2 (t)(5ij + hij)dx l dx^ , (2-18) 

the equations of motion for the dark matter density contrast d DM (k,t ) and velocity diver¬ 
gence #dm (k,t) for Fourier mode k are (e.g., [44]) 


^DM — 

$DM = —HO dm 


Oom h \ /3 ■ 

— + 2 + M^’ 


P (k 2 r , ; 

— 0(p — 0C7 DM 


M P1 \ a 


(2.19) 


and similarly for the scalar field perturbations, 


5(f) + 3 H5cj) + I — + Vrf (j> 


1 + 2 ^^ _ ~~ Pompom , 


( 2 . 20 ) 


where h = h l i is the trace of metric perturbation hjj. Note that in non-interacting models the 
dark matter velocity divergence 0 DM is exactly vanishing at all times because the synchronous 
coordinates are defined such that initially 0J5 m = = 0, and there is no source term to 

perturb (9 DM away from zero. In chameleon models, however, the DM-DE coupling engenders 
a source term proportional to the scalar field perturbation dcj), so that it becomes necessary 
to track the evolution of d DM as well. 

The equation of motion for the metric perturbation h, obtained from the perturbed 
Einstein equation, likewise needs to be modified according to 


i a 2 Hh = 

a A 

2 M|j 

II 

•c- 

CN 

a / 

2Mp 2 , ^ 

k 2 

- 2 —r) = 
a z 

3 1 

Ml 


pA + dp<t, I , 

yy (pi. + Pi)0i + ak~ 

yy 5Pi+~ v^dcj) j, 


Z=7,1/,B,DM 




( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


where dp^ = (f)6(f) + V^Scj), and we observe that perturbations in the cj) field do not contribute 
to the traceless part of the spatial stress tensor. Note that the evolution equations presented 
in this section are common to all scalar-field-based couple DM-DE models; the uniqueness of 
the chameleon model lies in the behaviour of the potential and the existence of an attractor 
solution. 
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It is instructive to rewrite equation (2.19) as a second order differential equation for <5 DM , 

^ + 2H^ + l(h + 2 (2.24) 

where we have ignored terms proportional to (f> because of the slow-roll conditions (2.17), and 
solved equation (2.20) for Scj> assuming a steady state (i.e., 6(j) = 6<j) = 0). While the term 
dependent on the metric perturbation h arises from standard Einsteinian gravity and can be 
easily constructed from equations (2.21) and (2.23), the additional term proportional to f3 2 
can be interpreted as a fifth force acting on the dark matter perturbations. The (comoving) 
range of this fifth force is determined by 

A P (a) = a" 1 ^ 2 , (2.25) 

so that only those scales satisfying k > A F 1 will feel the force’s effect. For potentials of the 
form (2.8), we find 


A F (a) = a 


-l 


cj) a+ 2 


a(a + 1 )M, 4 M P ° 


a 2 /7 ( 7 1 , 




a 


02 

“dm 


3(3 2 


= —/ 3 ( 0 — <fro)/Mpi 


(2.26) 


where we have assumed in the last equality a« 1. Thus, while A P today is comparable to 
the present-day Hubble length, at early times the range of the force is strongly suppressed 
by the scale factor a. At the time of CMB decoupling, for example, equation (2.26) evaluates 
to _ 


A F (a*) 


a* / 2 (a*ff *) _1 


a 


n T 


3/3 2 


(2.27) 


which shows that A P (a* ~ 0.001) is no more than about ICC 5 times the comoving Hubble 
length at a*, and consequently completely out of the observable range of the CMB primary 
anisotropies (> 0(0.1) x (a*i7*) -1 ). 


3 Effects on cosmological observables 

We describe in this section the observational consequences of the chameleon model for the 
current generation of cosmological probes. In particular, we discuss the parameter degenera¬ 
cies expected to arise between the a and (3 and other cosmological parameters, and determine 
which combinations of cosmological observations would be capable of lifting these degenera¬ 
cies. We defer the discussion of chameleon effects for future observations to section 5. 


3.1 CMB anisotropies 

Because the DM-DE coupling causes the dark matter density to dilute more slowly than a -3 , 
for the same present-day value, the chameleon model has a lower dark matter density at 
high redshifts than does the standard ACDM case. This has the effect of delaying the epoch 
of matter- radiation equality as per 


w DM e“ Wo/Mpl + w B 1 

w 7 1 + 0.2277V eff 

w D Me~“( ,i2 /“ DM ” 1 ) + uj B 1 

^ 1 + 0.227-/V e fT ’ 


(3.1) 
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where N eS = 3.046, and the second approximate equality follows from equation (2.15). 
Changing the epoch of equality has a profound impact on the odd acoustic peak ratios of the 
CMB temperature anisotropies, and with seven acoustic peaks now observed by Planck, z eq 
and hence cu DM e _Q ^ 2/,aJDM ^ can be considered a well-constrained quantity. 2 

However, merely measuring cu DM e _a ^ / WDM_1 ) is clearly not sufficient to determine all 
of the chameleon model parameters, because of the intricate dependence of the exponent on 
the Hubble parameter h, the potential’s slope a, as well as w DM itself, let alone the missing /?. 
Breaking the degeneracy requires that we measure three other parameter combinations; two 
of these come automatically from the CMB temperature anisotropies: 


Angular sound horizon. The 9 S = r*/D A parameter determines the positions of the 
CMB acoustic peaks. In the numerator, r* = r s (a*) = dtc s (t)/a(t) is the comoving sound 
horizon at the time of CMB decoupling. For the chameleon model this evaluates, up to a 
constant factor, to 


h a<iq In 

•\/l + R(z *) + \J R(z *) + R(z eq ) 

/ 3 (woMe-^/^DM-i) + u B )R(z eq ) “ 

1 + \/ ^(^eq) 


(3.2) 


where R(z) = (3/4)(w B /w 7 )a. In the denominator, D A = Z) A (a*) is the angular diameter 
distance to the last scattering surface, which in the chameleon model takes the approximate 
form 


D* 

IS A 


oc 


f 

J a* 


da 


1 2 \j (w DM + w B )a“ 3 + (h 2 - w DM - w B )e 


(3.3) 


up to a constant factor. Measuring z eq and u B from the peak ratios completely fixes r*, but 
not D a which depends on a different combination of cj dm and h, as well as a through the 
effective dark energy equation of state w eH (as already discussed in section 2.2, w efi does not 
depend on f3). Therefore, 9 S contributes towards breaking the (a, h, w DM )-degeneracy. 


Late Integrated Sachs Wolfe (ISW) effect. The late ISW effect is primarily sensitive 
to the evolution of metric perturbations on those length scales that enter the Hubble horizon 
while the universe transitions from matter to dark energy domination at low redshifts, z < 1. 
The metric perturbations evolve according to equations (2.21) to (2.23) in the chameleon 
model, and have contributions from both the matter component and the scalar field. 

The scalar field density perturbation bp^ is given approximately by 

<5/30 ~ V^bcj) « 3/? 2 fl DM (t) £2^ i _|_ a 2 V^jk 2 ^ DM ’ 


where we have used = —(/3/M Pl )p DM , and assumed the steady-state solution of equa¬ 
tion (2.20) for bcf> (i.e., 5d) = b(f> = 0). Feeding this solution into equations (2.21) and (2.23) 
also allows us to rewrite the evolution equation (2.24) for the matter perturbations in the 
low-redshift universe as 


<5dm + 2 H5 om — —H^ fI E 


3 

— J 
2 


i(t) 


1 + 


2/3 2 


1 + a^V H /k 2 


1 312dm (f) 


a 2 H 2 {t) 

k 2 


-'DM 5 


(3.5) 


where we have neglected the baryon component for simplicity. Since according to equa¬ 
tion (2.26) the range of the fifth force has recently attained Hubble length, we can establish 

2 The baryon density wb is independently well constrained by the odd-to-even acoustic peak ratios. 
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Figure 2. Bardeen potential 4>#(/c) as a function of the scale factor a for k = 0.0001 /iMpc^ 1 
(black), 0.01 /iMpc -1 (red), and 0.1 hMpc^ 1 (cyan). Left: Chameleon models with a fixed a = 0.1, 
and /? values of 0.01 (solid), 0.05 (dashed), 0.1 (dot-dashed) and 0.2 (dotted); the bottom portion 
shows the fractional difference A$#/$# between the /3 = 0.01 and 0.2 cases. Right: Chameleon 
models with a fixed /3 = 0.1, and a values of 0.01 (solid), 0.05 (dashed), 0.1 (dot-dashed) and 0.2 
(dotted); the bottom portion shows the fractional difference between the a = 0.01 and 0.2 cases. 
In both panels all other cosmological parameters have been fixed to the ACDM best-fit values from 
Planck-)-WP+HST. 


that aV^/fc 2 < 1 for all observable scales at z < 1. It then follows from equation (3.4) 
that bp^ is strongly dependent on the coupling parameter /3, and its impact on the metric 
perturbations—both directly and indirectly through <f DM —is greatest on scales close to the 
Hubble length due to its a 2 H 2 {t)/k 2 dependence. We illustrate this point in figure 2 with 
the gauge-invariant Bardeen potential, 


$h = 


2 M 2 k 2 


E 




3 aH 

Pibi H- \Pi + Pi)Vi 


+ bp<j> + 3 a~Hf>5(f> 


(3.6) 


as a function of the scale factor a at k = 0.0001,0.01,0.1 /i Mpc" 1 for various combinations 
of a and (3. 

A second interesting feature in equations (3.4) and (3.5) is that, while the low-redshift 
behaviour of H(t ) can be conveniently reparameterised in terms of an effective dark energy 
equation of state w eS , the time-dependent reduced density parameter 0 DM (t), given in this 
case by 




Pnha 3 


PmU 3 + P<t>e ^ M/M p1 


(3.7) 


does not share the same mapping, i.e., j n the demonimator does not equate 

to /3 DB of equation (2.12). Thus, even if the f3 2 term should be absent in equation (3.5), the 
evolution of <f DM and hence &h in the chameleon model cannot be mapped to that in a non¬ 
interacting scenario with a dark energy equation of state w eS . In other words, the evolution 
of d DM and &h in the low-redshift universe probes a combination of (a, h,uj DM ) different from 
that determined by w eS . 


Note that while measurements of the CMB angular damping scale 9 d = r* A /D* A by 
ACT [31] and SPT [32] have been instrumental in procuring the first CMB-only constraints 
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on the radiation energy density, for the chameleon model 9 d offers no new information beyond 
what can already be obtained from the acoustic peaks. To see this, we note that the diffusion 
damping scale r* = r d (a*) can be approximated as [40] 


(2ir) 


f 


da 


R 2 + (16/15) (1+R) 


o a 3 aTn e H l 6(1 + R 2 ) 


oc 


(2vr) 2 

r* da 

'R 2 + (16/15) (1 + R)' 

V^D M e- a (' l2 /^DM-i) + Ub .7 

o a^/ a + a eq aTn e 

6(1 + R 2 ) 


(3.8) 


where ut is the Thomson scattering cross-section, and n e is the free electron number density. 
Clearly, after z eq and u B have been fixed, r* d is likewise completely determined, while the an¬ 
gular diameter distance D* A is probed by the angular sound horizon scale 9 a . The information 
contained in the angular diffusion scale 9 d is therefore redundant. Indeed, as we shall see 
in section 4, adding ACT and SPT data to WMAP’s acoustic peak measurements does not 
improve the constraints on the chameleon model parameters. 


3.2 Non-CMB observables 

One more orthogonal parameter combination needs to be constrained in order to eliminate 
all degeneracies. This can come from any one of the following observables: 


1. The Hubble expansion parameter h can be fixed directly by measurements in our local 
neighbourhood. 


2. The current generation of baryon acoustic oscillations (BAO) measurements constrains 
the angular scale r s (z d )/D v (z BAO ), where z d ~ 1020 denotes the redshift of the baryon 
drag epoch, z BAO the effective redshift at which the BAO is observed, and 


D v (z) 


z(l + z) 2 



H{z) 


(3.9) 


is a distance scale subsuming both the angular diameter distance D A (z) and the Hubble 
expansion rate 


H{a) cc (o; DM + w B )a 3 + {h 2 - w DM - w B )e 3 -A ; . (3.10) 

The BAO peak shares a common physical origin with the CMB acoustic peaks. How¬ 
ever, because the BAO measurement takes place at a low redshift (.Zbao < 1), the 
parameter combination to which it is sensitive is vastly different from that probed by 
the high-redshift r*/D A of equations (3.2) and (3.3). 

3. An alternative low-redshift distance constraint comes from measurements of the lu¬ 
minosity distance D L (a) versus redshift using type la supernovae. In a flat spatial 
geometry, D h (a) = D A (a)/a 2 . 

4. The large-scale matter power spectrum extracted from galaxy clustering surveys yields 
several useful pieces of information. Firstly, the “turning point” depends on the co¬ 
moving wavenumber 

k eq = a eq H(a eq ) « 4.7 x 1CT 4 /i - 1 ^/(w DM e" a ( ,l2 /“ DM “ 1 ) + w B )(l + z eq ) /iMpc -1 , (3.11) 
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which, after fixing z eq and w B with the CMB, effectively depends only on the Hubble 
parameter h (recall that k is measured in units of /iMpc -1 ). Secondly, BAO features 
are also present in the present generation of power spectrum measurements (e.g., [45]), 
which offers the same handle on a low-redshift distance scale as discussed immediately 
above. Thirdly, the fifth force induced by the DM-DE coupling in principle enhances 
clustering in a scale-dependent fashion even on scales well within the Hubble horizon 
according to equation (3.5), thereby distorting the matter power spectrum. However, 
because the range of the force has only reached Hubble length at recent times, we 
expect this subhorizon spectral distortion to be confined mainly to the nonlinear scales 
(k > 0.2 ZiMpW 1 ), where the density perturbations have had more time to evolve under 
the influence of the fifth force [29] and which are beyond the scope of our perturbative 
analysis. Nonetheless, we note that in certain f(R) or interacting quintessence models 
it is possible to shift the clustering enhancement to the linear scales (e.g., [46, 47]). 

Other low-redshift matter distribution measurements such as the cluster mass function 
[48-50], cosmic shear [51-53], and the Lyman-ct forest [54, 55] are also potentially useful for 
constraining chameleon cosmology in that they provide a direct measurement of the density 
fluctuation amplitude which should be sensitive to the non-standard low-redshift evolution 
of the matter density perturbations induced by the DM-DE coupling. However, to make use 
of these measurements requires nonlinear modelling, which is beyond the scope of this work. 

4 Cosmological data analysis 

Following the discussions of section 3, we consider four generic types of cosmological measure¬ 
ments: temperature and polar polarisation power spectra of the CMB anisotropies, the BAO 
scale, the large-scale matter power spectrum, and direct measurements of the local Hubble 
expansion rate. The details of the actual data sets used are related below. On these data 
sets we perform a Bayesian statistical inference analysis using the publicly available Markov 
Chain Monte Carlo (MCMC) parameter estimation package CosmoMC [56] coupled to the 
CAMB [42] Boltzmann solver modified in accordance with section 2. With the exception of 
the local Hubble parameter measurement, the likelihood routines and the associated window 
functions are supplied by the experimental collaborations. 

Table 1 shows the fit parameters used in the analysis and their associated (flat) priors 
ranges. We adopt the Gelman-Rubin convergence criterion R — 1 < 0.02 when generating 
our multiple Markov chains, where R is the variance of chain means divided by the mean of 
chain variances. 

4.1 Data sets 

We split our analysis into two stages, the first centred on the WMAP 9-year measurements 
of the cosmic microwave background temperature and polarisation anisotropies [57], and 
the second on the more recent temperature measurements from the first data release of the 
Planck mission [5]. We describe these and other auxiliary data sets below. 

CMB anisotropies. In the stage 1, we use the temperature (TT), E-polarisation (EE), and 
E-polarisation (BB) autocorrelation power spectra, as well as the temperature-E-polarisation 
(TE) cross-correlation power spectra from WMAP [5]. To this we add the TT power spectra 
from ACT and SPT [58-62], which cover the range 600 < i < 3000 (SPT) and 500 < 
l < 10000 (ACT). The ACT/SPT data include in addition the deflection power spectra 
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Parameter 

Description 

Prior 

ACDM 

OJB 

Baryon density 

0.005 -)■ 0.04 

- 

£Vdm 

Cold dark matter density 

0.01 0.99 

- 

h 

Hubble parameter 

0.4 —> 1.0 

- 

Zre 

Redshift of reionization 

3 35 

- 

n s 

Spectral index of the primordial power spectrum 

0.5 -> 1.5 

- 

ln(10 10 A s ) 

Amplitude of the primordial power spectrum at k = 0.05 Mpc -1 

2.7 4.0 

- 

a 

Slope of the chameleon potential 

10“° ->■ 5 

0 

/3 

Strength of the chameleon coupling 

10“ 5 -)• 5 

0 

Neff 

Number of effective massless neutrino families 

1 ->• 9.8 

3.046 


Table 1 . Fitting parameters in our MCMC analysis and their associated flat prior ranges. The six 
parameters in the first block are the standard ACDM parameters. The chameleon model parameter 
space subsumes the ACDM parameter space, and has an additional two parameters a and /3 which, 
when set to zero, reduce the chameleon model to ACDM. The last entry N eS parameterises the 
radiation excess; in ACDM this is fixed at 3.046. 


(4>(j>) due to gravitational lensing, which we also analyse in our study. In stage 2, following 
the guidelines of [5], we use the TT spectra from the Planck mission in conjunction with the 
WMAP 9-year polarisation measurements (WP). To this we add the ACT/SPT TT spectrum 
at high multipoles (HighL), which helps to constrain the foreground nuisance parameters used 
in the Planck data analysis. 

Local Hubble parameter measurements. We adopt the values of Hq measured in our 
local neighbourhood by the Hubble Space Telescope. In stage 1, we use Hq = 74.2 ± 3.6 km 
s^ 1 Mpc -1 [63], published in 2009, in conjunction with our analysis of the WMAP 9-year 
CMB measurements. This value of Hq was also adopted by the WMAP collaboration in their 
9-year analysis [57], and we opt to keep it to facilitate comparison. In stage 2, we adopt a 
more recent measured value, Hq = 73.8 ± 2.4 km s -1 Mpc" 1 [64], in combination with our 
Planck analysis. Because each value is used exclusively with either WMAP or Planck, we 
adopt the same acronym HST for both measurements. 

Large-scale matter power spectrum We use the red luminous galaxy survey from the 
Sloan Digital Sky Survey DR-7 (SDSS) [45] to constrain the matter power spectrum in the 
/c-range [0.02,0.2] /iMpc -1 , in conjunction with the WMAP 9-year data in stage 1 of the 
analysis. The power spectrum measurement contains both broadband information on the 
scale-dependence, as well as geometric information in the form wiggles from the baryonic 
acoustic oscillations. In combination with the WMAP CMB measurements, the latter infor¬ 
mation generally suffices to constrain model parameters that affect the late-time expansion of 
the universe, while the former is primarily useful for extended models in which the perturba¬ 
tions pick up an additional scale-dependence at late times (e.g., models with finite neutrino 
masses) [65]. In the context of chameleon models, we expect geometric information to be the 
more useful of the two (the additional scale-dependence in chameleon models are in any case 
all confined to the nonlinear scales beyond k ~ 0.2 AMpc -1 , as discussed in section 3). We 
therefore do not use the matter power spectrum in stage 2 of the analysis in combination 
with Planck data, and opt instead for direct measurements of the BAO scale. See next. 

BAO scale. We use measurements of the BAO scale by SDSS-DR [66], SDSS-DR9 [67], 
and 6dFGS [68]. For cosmological models with no non-trivial scale dependence, geometric 
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Model 

Data 

a 

P 

Aeff 

WDM 

h 

AXa/ Adof 

Axb /Adof 

Chameleon 

W+HST 

< 0.23 

< 0.19 

3.046 

0 . 1098:099 

o.7o8:8l 

+ 0 . 0/2 

+ 0 . 1/1 

W+HST+ACT/SPT 

< 0.24 

< 0.19 

3.046 

0.1098:188 

0 . 708:88 

+ 1 . 3/2 

+0.7/1 

W+SDSS 

< 0.48 

< 0.19 

3.046 

0 . 1128 T 03 

0.658:88 

- 0 . 1/2 

-0.3/1 

W+HST+SDSS 

< 0.27 

< 0.19 

3.046 

o.ii 28 :J 88 

0 . 688:88 

+ 0 . 1/2 

+ 1 . 1/1 

W+SNIa+SDSS 

< 0.19 

< 0.19 

3.046 

o.ii38:188 

0 . 688:88 

- 0 . 4/2 

-0.3/1 

P+WP+BAO 

< 0.17 

< 0.19 

3.046 

0.1188:118 

0.678:88 

+0.7/2 

+0.5/1 

P+WP+HST 

< 0.09 

< 0.19 

3.046 

0.1178:118 

0 . 688:88 

+0.4/2 

+4.7/1 

P+WP+HighL+BAO 

< 0.17 

< 0.19 

3.046 

o.ii88:118 

0.678:88 

-0.7/2 

+0.4/1 

P+WP+HighL+HST 

< 0.09 

< 0.19 

3.046 

o.ii78:118 

0.678:88 

- 0 . 6/2 

+4.7/1 

P+WP+HighL+BAO+HST 

< 0.08 

< 0.19 

3.046 

o.ii78:lli 

0.678:88 

- 0 . 6/2 

+3.1/1 

2 

g 

< 

a? 

W+HST 

- 

- 

A q 6.1 

4 -<J 2.8 

o.i38:18 

0 . 758:88 

- 0 . 1/1 


W+HST+ACT/SPT 

- 

- 

q r)3.9 

d.^ 2.6 

o.i48:ll 

0 . 748:88 

+ 0 . 6/1 


W+SDSS 

- 

- 

4.55:? 

o.i48:18 

0 . 758:88 

+0.2/1 


W+HST+SDSS 

- 

- 


o.i48:l? 

0 . 748:88 

-1.0/1 


W+SNIa+SDSS 

- 

- 

4.65:1 

o.i48:?8 

0.768:88 

-0.1/1 


P+WP+BAO 

- 

- 

q q3.9 

d.02.8 

0.1248:111 

0.698:88 

+0.2/1 


P+WP+HST 

- 

- 

3.71:5 

0.1288:118 

0 . 738:88 

-4.3/1 


P+WP+HighL+BAO 

- 

- 

q q3.8 

0 . 0 2.7 

0.1238:118 

0.658:88 

-1.1/1 


P+WP+HighL+HST 

- 

- 

q r)3.ti 

O.Z 2.7 

o.ii98:l?? 

0 . 688:88 

-5.3/1 


P+WP+HighL+BAO+HST 

- 

- 

q 

O.Z 2.7 

0.1208:118 

0 . 688:88 

-3.8/1 



Table 2. ID marginal statistics for the chameleon parameters a and fj in the chameleon fit (top 
block) and for N eB in the lV eff ACDM fit (bottom block) derived from various data combinations. The 
abbreviation “W” stands for WMAP, while “P” is Planck. For one-sided intervals we show only the 
95% credible limit, while for two-sided intervals we show the 95% central credible interval (see [70] for 
definition) as well as the posterior mean. The second last column compares the model’s best-fit x 2 
value relative to the ACDM best-fit as per equation (4.1) in the face of A dof additional degrees of 
freedom, while the last column contrasts the chameleon and the -/V eff ACDM fits as per equation (4.2). 


information extracted from large-scale structure surveys suffices to constrain the model pa¬ 
rameters [65], and has the added advantage of being much less prone to the nonlinearity 
issues that plague broadband measurements of the large-scale matter power spectrum. 

Type la supernovae. We include the supernovae data set Union2 compilation [69], la¬ 
belled “SNIa”, in stage 1 of the analysis. The standardised luminosity of Type la supernovae 
is a classic measure of the redshift-dependence of the (h-normalised) luminosity distance. 
This compilation consists of 557 supernovae collected from different surveys, and was a stan¬ 
dard compilation used also by the WMAP collaboration in their 9-year data analysis [57]. 

4.2 Parameter constraints 

Table 2 summarises the constraints on the chameleon model parameters a and /3 from 
our MCMC analysis using various data combinations. The corresponding two-dimensional 
marginal posteriors in the (a, /3)-plane are displayed in figure 3. 

Clearly, none of the data combinations prefer the chameleon model over the basic ACDM 
model. This can be discerned firstly from the fact that only one-sided (instead of two-sided) 
limits exist for the chameleon parameters a and /3, indicating that a = (3 = 0 are completely 
compatible with all data combinations. 3 Secondly, the different between the best-fit x 2 values 

3 Note that in practice, because we use equations (2.15) and (2.16) to determine the present-day field 
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Figure 3. 2D marginal 68% and 95% contours in the (a,/3)-plane derived from various data com¬ 
binations. Left: WMAP temperature and polarisation measurements combined with HST (solid 
black), SDSS (dashed red), HST+SDSS (dotted blue), HST+ACT/SPT (dot-dashed magenta), and 
SDSS+SNIa (solid cyan). Right: Planck temperature and WMAP polarisation data combined with 
HST (solid black), BAO (dashed red), HighL+HST (dotted blue), HighL+BAO (dot-dashed ma¬ 
genta), and HighL+HST+BAO (solid cyan). 


of the chameleon and the ACDM model, 


^Xa — Xmodel XACDM) (^-l) 

is at best —0.66 for the Planck+WP+HighL+BAO data combination. This is barely beyond 
the reliability limit of |Ay 2 | > 0.6 of the BOBYQA bounded minimisation routine [71] used 
to find the best-fit y 2 values (the Metropolis- Hastings algorithm in standard MCMC analyses 
is not suitable for finding the maximum of a multi-dimensional likelihood function), and is 
in any case too small an improvement to the fit to warrant the introduction of two extra fit 
parameters. 

In terms of the limits on the chameleon parameters a and (3, there is virtually no dif¬ 
ference between data combinations in their constraining power on the coupling parameter (3; 
all return the same 95% upper limit (3 < 0.19. This can be understood as follows: as already 
discussed in sections 2.2 and 3.1, neither the effective dark energy equation of state w eff nor 
the high-redshift observable z eq is sensitive to /3; the effect of /3 is felt only through the late 
ISW effect, because of the enhanced Bardeen potential &h at 2 < 1 on close to the Hubble 
length as shown in figure 2. Figure 4 shows the CMB TT spectrum of various chameleon 
models relative to the ACDM best-fit to the Planck+WP+HST data combination. Although 
the /3-dependence is clearly strong, the late ISW effect evidently manifests itself only at low 
multipoles i < 30. This means firstly that, besides the CMB TT spectrum, the cosmological 
observations considered in this analysis play no role in constraining j3. Secondly, swapping 
WMAP for Planck CMB temperature data, both of which are limited only by cosmic variance 

value <f >o and the field mass M^, it is necessary to impose a nonzero lower limit on the priors on a and /?, both 
chosen here to be ICC 5 as shown in table 1, in order to avoid an artificial divergence at a = f) = 0. However, 
because the phenomenology of a = /3 = 10 -5 is for all purposes indistinguishable from ACDM, we choose not 
to enforce the distinction between the two cases. 
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Figure 4. CMB TT spectrum of various chameleon models relative to the ACDM case. The ACDM 
spectrum corresponds to the best-fit model of the data set Planck+WP+HST, while the chameleon 
models have the common parameters (a,L 0 VM , h) = (0.1,0.117,0.688) and, from bottom to top, /3 = 
0.05,0.1,0.15,0.2,0.25,0.3. 

at the low multipoles, also cannot improve the sensitivity to fi. Further improvements can 
come about by either looking at the effects of the fifth force induced by /3 on very nonlin¬ 
ear scales, or, in the future, the impact of the /3-dependent metric perturbations on full-sky 
cosmic shear surveys such as the ESA Euclid mission [72] 

For the slope parameter a, the constraint generally improves when WMAP data are 
swapped for Planck measurements. This can be understood as a consequence both of a more 
well-constrained z eq and of a more accurate determination of 0 S . The latter in particular 
constrains the effective equation of state w eS , whose departure from the canonical value of — 1 
is now, as shown in figure 5 and table 3, restricted to no more than 20% (at 95% confidence) 
at 2 < 0.5 and 50% at z < 0.8. Adding the HST measurement of Hq to Planck improves the 
constraint on a even further (almost a factor of two). However, this improvement may be an 
artefact originating from the 2a+ discrepancy between the HST measurement and the h value 
inferred from Planck data in the minimal ACDM model [5, 6]. Indeed, keeping u; DM fixed, we 
see, on the one hand, from equation (3.1) that a must decrease with an increasing h in order 
to maintain the same z eq . On the other hand, from equation (3.3), maintaining the same D* K 
in the face of a smaller a requires that we reduce h. Thus, the net effect is that a goes down, 
while the inferred value of h changes only marginally, as suggested by the numbers in table 2. 
This incompatibility between data sets is also confirmed by the 2D marginal posterior in the 
(a, i7o)-plane in figure 6, which shows a clear preference for a negative a value—which is 
not within the chameleon parameter space—if Hq was indeed as large as 73.8 km s -1 Mpc~ x 
preferred by HST. 

Comparing our results with previous investigations of interacting DM-DE models, our 
limits differ in significant ways. The analyses of [29, 73] dealt with a scenario identical to 
ours, and it was argued in these works that /3 ~ 0(1) could be compatible with observations. 
However, this conclusion was based on the argument that with fi and a fixed at 1 and 0.2 
respectively, a best-fit y 2 “close to” the ACDM best-fit could be obtained when the usual 
ACDM parameters are let to vary freely. Our analysis here, in contrast, allows all parameters 
to vary, so as to locate the values of a and f3 actually preferred by the observational data. 


\ 

\ 

\ \ 


\ \ 

N \ \ 
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Planck+WP+HST 


Planck+WP 





-1.1437 -1.1078 -1.0719 -1.036 -1.0001 

W eff(°‘2) 



0 0.05 0.1 0.15 0.2 0.25 0.3 



-1.2094 -1.1571 -1.1047 -1.0524 

"JO-3) 


-1.2455 -1.1842 -1.1228 -1.0615 -1.0002 

w JO-3) 




-2.3815 -2.0362 -1.6908 -1.3455 -1.0001 

w J°- 8 > 



-2.7664 -2.325 -1.8836 -1.4422 -1.0008 

w JO-8) 


Figure 5. 2D marginal 68% and 95% contours in the (a, /3)-plane derived from Planck+WP+HST 
(left) and Planck+WP (right), and the corresponding effective equation of state for the dark energy 
(coloured dots) at four different redshifts, from top to bottom, z = 0.1,0.2,0.3,0.8. 


References [74] and [75] also considered a coupled scenario in which the scalar potential V((f)) 
is of the runaway form (2.8), and obtained constraints on (5 that are considerably tighter 
than ours. We note however that the coupling function adopted in these analyses is of the 
form f(4>) = which differs from ours by a crucial minus sign, signifying that their 

effective potential V eS (cf>) does not possess a local minimum, nor < ■pit') an adiabatic attractor 
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z 

Weff(z) 

Planck+WP+HST 

Planck+WP 

0.0 

> -1.003 

> -1.005 

0.1 

> -1.03 

> -1.06 

0.2 

-1.06 

> -1.12 

0.3 

-1.11 

> -1.21 

0.4 

-1.16 

-1.32 

0.5 

> -1.2 

> -1.5 

0.6 

> -1.3 

> -1.7 

0.7 

> -1.4 

> -2.0 

0.8 

> -1.5 

-2.4 


Table 3. ID marginal 95% lower limit for the effective dark energy equation of state u> e ff in the 
chameleon model at various redshifts, derived from the data combinations Planck+WP+HST and 
Planck+WP. 


solution. For this reason, the results of [74] and [75] cannot be meaningfully compared with 
ours in the context of parameter inference. 

4.3 Chameleon mimicking dark radiation? 

As discussed at the beginning of section 3, the nonstandard time-evolution induced for the 
dark matter density by the DM-DE coupling tends to delay the epoch of matter-radiation 
equality relative to the standard ACDM case with the same w DM . This raises the possibility 
that chameleon phenomenology might mimic the effect of dark radiation, which also has 
shifting z eS as its primary effect. We test this possibility by fitting the ACDM model extended 
with a free N eS to the same data combinations explored in section 4.2 for the chameleon model, 
and compute in each case the constraints on N eS as well as the best-fit y 2 value relative to 
ACDM as defined in equation (4.1). In addition, for each data combination we evaluate the 
y 2 difference between the best-fit chameleon model and its lV eff ACDM counterpart, defined 
as 

^chameleon ^ AT e ff ACDM * (4.2) 

This Ax 2 value will tell us whether or not chameleon models are able to mimic dark radiation. 

Table 2 shows the results of this exercise. Note that in deriving these constraints we have 
made use of the BBN consistency relation [76] in CosmoMC, which automatically adjusts 
the helium mass fraction Y He4 according to big bang nucleosynthesis predictions for each set 
of input parameters N eff and u B . Clearly, only the combinations WMAP+HST+SDSS and 
Planck+WP+HST show a 2cr+ preference for N eS > 3.046, and in the latter case, this is 
accompanied by a dramatic increase in the inferred h value ( h = 0.73(jgg) relative to other 
Planck data combinations. In terms of the best-fit y 2 , we find Ay 2 values of —0.98 and —4.32 
respectively relative to ACDM, suggesting a mild preference for a radiation excess especially 
in the latter case. 

Comparing the chameleon model with iV eff ACDM for these same two data combinations, 
we find the chameleon scenario consistently the worse performer of the two. The best- 
fit y 2 differences between the two models as per definition (4.2) are +1.13 and +4.68 for 
WMAP+HST+SDSS and Planck+WP+HST respectively. This can be traced to the fact that 
phenomenologically, the a parameter needs to be negative in order to accommodate the large 
Hq value preferred by HST (see figure 6 and discussion in section 4.2). A negative a, however, 
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Figure 6. 2D marginal posterior in the (a, iJo)-plane derived from Planck+WP+HST. 

does not belong to the chameleon parameter space. Thus, despite a passing similarity at the 
equality epoch, the degeneracy between Hq and the chameleon model parameter a ultimately 
works in the wrong direction to be able to resolve the tension between the HST measurement 
and the Planck-inferred Hq value. 

At this point, we note that our chameleon model is distinct from early dark energy 
models which do mimic a radiation excess because the dark energy equation of state in these 
models evolves to ~ 1/3 during radiation domination [77]. In contrast, the attractor solution 
in chameleon models ensures that the scalar field remains slowly rolling while sitting in 
the effective minimum of the potential, and therefore maintains a cosmological constant-like 
behaviour at all times. 

5 Future tests 

The two main predictions of the chameleon model are the apparent phantom behaviour of 
the effective dark energy equation of state w eS and a fifth force that enhances the clustering 
of dark matter as well as the gravitational potentials. The former is, for any nonzero /3, solely 
a function of the parameter a , while the latter is strongly dependent on j3. 

As we have seen in figure 5 and table 3, the effective dark energy equation of state is 
already constrained at 95% confidence to w eff > —1.2 at z = 0.5 and < —1.5 at z = 0.8 by 
Planck+WP+HST. It would therefore appear at first glance that the Dark Energy Survey 
(DES), with la sensitivities a(w o) ~ 0.078 and a(w a ) ~ 0.3 [78] for the parameterisation 

w(a) = w 0 + w a (l - a), (5.1) 

may not be sufficient to further tighten the constraint on the chameleon parameter a. The 
ESA Euclid mission on the other hand, with a(wo) ~ 0.018 and cr(w a ) ~ 0.15 [72], should 
produce some improvements. 

We caution however that this interpretation is not strictly correct: parameter sensitiv¬ 
ities derived for a simple dynamical dark model such as equation (5.1) do not automatically 
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translate to bounds on the chameleon model parameters. This is because ongoing and future 
large-scale structure surveys will derive their constraints on the dark energy equation of state 
from both distance measures and the growth of density perturbations; apart from its effect 
on the late ISW effect, the latter observable has not been used in the present analysis of 
existing cosmological data. Recall from our discussions in section 3.1 that in the chameleon 
scenario, the evolution of the dark matter density parameter H DM (f), equation (3.7), which 
governs the gravitational potential felt by the density perturbations, maps to an effective 
dark energy equation of state different from the canonical w eS responsible for the background 
Hubble expansion. This additional feature should allow us to distinguish cleanly between 
simple dynamical dark energy and interacting/chameleon models. A detailed analysis of how 
well future surveys will be able to select between dark energy models is beyond the scope 
of the present work, but we note that publicly available forecast codes such as [79-81] exist 
which, with minor modifications, will be equal to the task. 

As for the fifth force, there are two possible places to look for its effects. One is the 
very nonlinear scales, where the perturbations have had the most time to evolve under the 
force’s influence. Given that the ft parameter is already constrained by the late ISW effect 
to less than 0.2, we expect the fifth force to manifest itself on subhorizon scales, according 
to equation (3.5), as a maximum 8% correction to the (scale-dependent) effective Newtonian 
constant; nonlinear evolution should further enhance the effect (although difficult to predict 
precisely). The second possibility is to take advantage of the almost full sky coverage of future 
surveys such as the ESA Euclid mission. In the same way that the scale-dependent evolution 
of the Bardeen potential at z < 1 and k ~ ( aH ) _1 induced by [3 (see figure 2) enhances 
the late ISW effect, we expect an analogous distortion of, e.g., the cosmic shear convergence 
angular power spectrum, at similarly low multipoles l < 30. Of course, detection of this 
signature and ultimately the extent to which one could improve upon existing constraints 
on ft will still be, like the late ISW effect, limited by cosmic variance. To address this question 
quantitatively we would need to perform a parameter sensitivity forecast, which we defer to 
a future work. 

6 Conclusions 

We have considered in this work a simple interacting DM-DE/chameleon model in which 
a scalar field slowly rolls and adiabatically tracks the minimum of an effective potential 
starting from deep in the radiation-dominated era. We have tested this model against recent 
cosmological observations, especially new CMB temperature measurements from the Planck 
mission, and have found constraints on the two model parameters: a, slope of the scalar 
potential, and /3, the coupling strength. We find that while the constraint on the coupling 
strength ft has remained unchanged between WMAP and Planck, essentially because this 
parameter is constrained primarily by the late ISW effect, the upper limit on a, which has 
the dual effect of shifting the matter-radiation equality epoch and altering the effective dark 
energy equation of state, has improved somewhat. 

Because of its potential to shift the epoch of matter-radiation equality, we have tested 
also the possibility that such a chameleon model might mimic the phenomenology of dark 
radiation, and resolve the apparent tension between the relatively large Hubble expansion 
rate measured in our local neighbourhood by HST and the low value inferred from the Planck 
CMB data. We find however that the degeneracy between a. and Hq in the chameleon model 
works in the wrong direction: a larger Hq value in tends to push a towards zero, i.e., the 
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cosmological constant limit. In fact, the ultimate effect of the HST/Planck tension on the 
chameleon model is to put an artificially tight constraint on a. 

Taking into account the allowed model parameter space, we have discussed possible 
signatures of this type of chameleon model for ongoing and future large-scale structure surveys 
such as DES [78] and the ESA Euclid mission [72], The salient effects are (i) a mismatch 
between the effective dark energy equation of state inferred from the background expansion 
and from the growth function which should allow us to determine or further constrain the 
value of a, and (ii) distortion of the large-scale distribution power spectrum on very nonlinear 
and close-to-horizon scales, both of which are predominantly sensitive to (3. To determine 
the precise sensitivities of DES or Euclid to the chameleon model parameters would however 
require a dedicated forecast analysis. We leave this exercise for a future work. 
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